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Abstract 

The procedure for obtaining integrable vertex models via reflection matrices on the 
square lattice with open boundaries is reviewed and explicitly carried out for a number 
of two- and three-state vertex models. These models include the six-vertex model, the 
15-vertex A^ model and the 19-vertex models of Izergin-Korepin and Zamolodchikov- 
Fateev. In each case the eigenspectra is determined by application of either the algebraic 
or the analytic Bethe ansatz with inhomogeneities. With suitable choices of reflection 
matrices, these vertex models can be associated with integrable loop models on the 
same lattice. In general, the required choices do not coincide with those which lead to 
quantum group-invariant spin chains. The exact solution of the integrable loop models 
- including an 0(n) model on the square lattice with open boundaries - is of relevance 
to the surface critical behaviour of two-dimensional polymers. 
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1 Introduction 



Vertex models are usually solved with periodic boundary conditions |[ [5[ |]]. Pe- 
riodic boundaries ensure that the transfer matrix r(u), formed as the trace over the 
auxiliary space of the monodromy matrix T(u), forms a commuting family if T(u) is 
constructed from solutions R(u) of the Yang-Baxter equation. The consequences of 
this commutativity are well-known, and allows exact results to be obtained for various 
physical quantities related to the models. 

The coordinate Bethe ansatz -for vertex models- (see e.g. Ref. 0), of course, pre- 
dates Baxter's commuting transfer matrices. While simple enough for two-state models 
like the six-vertex model, it becomes too cumbersome for models with a higher number 
of states allowed on the edges. This is to be contrasted with the quantum inverse 
scattering method || and related algebraic and analytic Bethe ansatz techniques, which 
rely on the commutativity of t(u) in order to diagonalize it. Nevertheless, the first case 
of a vertex model with non-periodic boundaries was solved using the coordinate Bethe 
ansatz ||. This was for the six- vertex model on the lattice C, being the square lattice 
rotated by 45° and with open boundaries on two sides (Fig. 1).[] Subsequently |J 
this has been extended to solve a three-state vertex model on the same lattice. The 
fact remains, however, that the method is somewhat unsystematic, especially for open 
boundaries. 

On the other hand, Sklyanin J7| has given a construction of commuting transfer 
matrices t(u, a?), out of which he obtained the Hamiltonian for the spin-i XXZ chain 
with open boundaries (which has also been solved using the coordinate Bethe ansatz 
[[|, ||). This transfer matrix t(u, 10) was constructed in the spirit of the quantum 
inverse scattering method and involves the use of reflection or K-matrices which satisfy 
a boundary version of the Yang-Baxter equation, first discovered in the framework of 
factorized S'-matrices || pTOfl . However t(u,u>) seemed a little mysterious as far as a 
vertex model interpretation was concerned, and it was not until fairly recently that 
it was shown |TT| that by specializing the inhomogeneities u>i in a specific manner it 
becomes precisely the transfer matrix toiu) for a vertex model on £ (i.e. the "diagonal- 
to-diagonal transfer matrix" with open boundaries; see Fig. 2). Therefore the Sklyanin 
scheme "explains" why the coordinate Bethe ansatz works for open boundaries in the 
same way that Baxter's commuting transfer matrices t(u) "explains" why it works for 
periodic boundaries. 

By now the Sklyanin scheme has been generalized |T2, [13], [14| to handle most of the 
i?-matrices constructed from (affine) non-exceptional Lie algebras [IB, [T(J. For some 
cases [[T7|, [13], |19[ the transfer matrix t(u,oj) has been diagonalized with modified 

1 We will sometimes refer to this lattice as the square lattice with open boundaries. 
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Figure 1: The lattice £, of "extent" M and N (M = 10 and N = 6 in the fig- 
ure drawn) in the vertical and horizontal directions, respectively. Boundary 
conditions are open in the horizontal direction and periodic in the vertical. 



versions of the algebraic or analytic Bethe ansatz, the emphasis always being on the 
corresponding quantum chain with open boundaries. In this paper we use the Sklyanin 
scheme to obtain the corresponding integrable vertex models on the lattice L and their 
Bethe ansatz solutions. We will concentrate on several two- and three-state models, 
both for simplicity and also because these appear to be the most interesting physically. 
In particular we have in mind an application of these solutions to loop models [20] on 
the lattice C, especially those which have a bearing on the configurational statistics of 
polymers. 

To be precise, the vertex models concerned will turn out to have bulk weights Wi, 
i = 1, ... ,19 given in Fig. 3 and boundary weights wf and wf, % = 1, 2, 3 given in Fig. 
4.0 The vertex model is defined through the partition function 

^vertex = E II (Boltzmann weights) , (1.1) 

configs vertices 

in the usual way. With periodic boundary conditions in the vertical direction, this 
partition function is equivalent to the trace of a product of transfer matrices tu(u), 
which then becomes the central object to study. 

The outline of the paper is as follows: In Section 2 we explain Sklyanin's method 
in greater detail and show how to obtain integrable vertex models on the lattice C In 



For a two-state model, Wi, i = 1, . . . , 13, and are necessarily zero. 
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Figure 2: The "diagonal-to-diagonal" transfer matrix tjj{u) which takes the 
state (ai, . . . , o^) to the state . . . , (3^)- All internal edges are summed 
over. 

the following few sections we carry out the procedure explicitly and obtain integrable 
open boundary versions of the six- vertex model, the Zamolodchikov-Fateev 19- vertex 
model, the model and the Izergin-Korepin model, respectively. The corresponding 
Bethe ansatz solutions are also given. In Section 7 we show how some of these models 
on L have an interpretation as loop models on the same lattice. Bethe ansatz solutions 
for these integrable loop models with open boundaries are given together with their 
intended applications. Finally we conclude in Section 8 with a discussion of several 
open problems related to vertex models on C. 

2 i^-matrices and the diagonal- to- diagonal transfer 
matrix 

Let us first review the case of integrability in the bulk, i.e. for vertex models on a 
standard square lattice with periodic boundary conditions. Such an integrable n-state 
vertex model is specified by an i?-matrix R(u) which acts on (D n ® C n and satisfies the 
Yang-Baxter equation 

Ri 2 (u)R 13 (u + v)R 23 (v) = R 23 {v)R l3 {u + v)R 12 (u), (2.1) 

where R 12 (u), R 13 (u) and R 23 (u) act on (D n <g)C"<g)C n , with R 12 (u) = R(u)®l, R 23 {u) = 
1 ® R(u) etc. The vertex weights are determined from the matrix elements of R(u) as 
shown in Fig. 5. For a two-state vertex model, the states 1 and 2 will be associated 
with an up-arrow and a down-arrow, respectively, while for a three-state model, the 
states 1,2 and 3 will be associated with an up-arrow, an empty edge and a down-arrow, 
respectively. 
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Figure 3: Allowed bulk vertices for integrable two- and three-state vertex 
models studied in this paper. Vertex i has an associated Boltzmann weight 

Wi. 
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Figure 4: Allowed boundary vertices for the vertex models concerned. The 
left (respectively, right) boundary vertices have Boltzmann weights wf (re- 
spectively, wf). 

Define the monodromy matrices 

T(u,u) = R al {u + u) l )R a2 {u + u 2 )---R aN {u + u N ), (2.2) 
f(u,u>) = R Na (u- uj n )- ■ ■ R 2a {u- uj 2 )R\a{u- uji), (2.3) 

where a denotes the auxiliary space and the quantities u>i are "inhomogeneities" . Graph- 
ically, a labels the horizontal direction whereas 1, . . . , N label the columns of the square 
lattice (from left to right). Due to the Yang-Baxter equation ( |2.1| ), the monodromy ma- 
trix T(u, u)) satisfies the intertwining relation 

R\2{u — v)T (u, a?) T (v, u>) =T iy, u) T (u, u))R X2 {u - v), (2.4) 

1 2 

where T (u, <jS) = T(u, a?) <S> 1 and T (u, uj) — 1 ® T[u, uj) etc. It then follows that the 
row-to-row transfer matrix t(u, uj) = tr T(u,<jj), where the trace is over the auxiliary 
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Figure 5: Graphical depiction of a matrix element of R(u). The indices run 
from 1 to n for an n-state model. 

space, forms a commuting family: 

[t(u,uj),t(v,u>)] = 0. (2.5) 

To determine an integrable vertex model on the lattice C, it is first necessary that 
the bulk vertex weights be specified by an i?-matrix R(u) as described above - together 
with a +45° rotation (cf. Fig. 3). The boundary weights then follow, as we shall see, 
from .fT-matrices which satisfy boundary versions of the Yang-Baxter equation |7], [11 
We will require |TB[ that R(u) satisfies the following properties 



regularity 
unitarity 
PT — symmetry 



J R 12 (0) = p(0) 1/2 P 12 , (2.6) 
R 12 (u)R^(-u)=p(u), (2.7) 
V 12 R 12 (u)V 12 = R 12 (u) ht \ (2.8) 



crossm 



ing - symmetry : Ri 2 (u) =V R%(-u - rj) V 1 . (2.9) 



Here, p(u) is an even scalar function, V is the exchange operator defined by V{f ®g) = 
(9 ® /)> U denotes transposition in the space i, rj is the crossing parameter and V 
determines the crossing matrix M = V l V = M l . The crossing-symmetry can be 
replaced by a weaker symmetry ||14|| , as will be necessary when we consider the 
model in Section 5. We note that unitarity and crossing-symmetry together imply the 
useful relation 

M R\ 2 2 (-u - r]) M R%{u -rj) = p(u). (2.10) 

Furthermore, we see that unitarity implies T(u, u)T(— u, us) = Ilili p( u ~ ^i)- There- 
fore, up to a scalar factor, T(—u, uj) is the inverse of T(u, uj). Both monodromy matrices 
(|2.2|) and ( p.3|) are necessary in order to construct integrable models on the lattice C 
The boundary versions of the Yang-Baxter equation are given by 0, [L3[ 



1 2 
R 12 (u - v) K~ {u)R^ 2 t2 {u + v) K~ (v) 
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K~ (v)R 12 (u + v) K~ (u)R^ 2 {u - v), (2.11) 

11 i 2 

R 12 (-u + v) (K+Y 1 (u) M~ l R$ 2 (-u -v-2rj) M(K + f 2 (v) = 

(K + ) t2 (v) M Rn(-u -v-2rj) M~\K + )^ {u)R$ 2 {-u + i>),(2.12) 



1 2 

where K~ (u) = K~(u)®l, K~ (u) = \ ®K~{u) etc. In practice, if K~(u) is a solution 
of (pHQ then 

K + (u) = K-(-u-r}fM (2.13) 



is a solution of (|2.12|) . Define the "doubled monodromy matrix" U(u,iv) by 

U(u, w) = T(m, w)/r (u)T(u, u>), (2.14) 



with the multiplication in the auxiliary space a. By (|2.11|) we have the "intertwining 
relation" 

Rn(u -v)U (u, u)R^(u + v)U (v, w) = 

?7 (v,u>)R 12 (u + v) U (u,u:)R^ 2 (u-v). (2.15) 



The main result is the following: If the boundary equations ( 2.1 1|) and ( |2.12| ) are 
satisfied, then Sklyanin's transfer matrix ]7|, |I^] 

t(u, u>) = tr K + (u)U(u, w), (2.16) 

where the trace is again in the auxiliary space a, forms a commuting family: 

[t(u,w),t(y,u)] = 0. (2.17) 

The proof of commutativity |7j involves showing that the object 

tT 12 B 12 {u,v) ht2 A 12 {u,v), (2.18) 



where A 12 (u,v) and B 12 (u,v) are the left (respectively, right) hand sides of (|2.12|) and 



( |2.15| ), can be manipulated using unitarity (|2.7|), equation fl2.10| ), and the properties of 



the trace tr and transposition t« to obtain t(u, uj)t(v, u>) (respectively, t{v , <-j)t(-u, uj)) 
times the scalar factor p(u—v)p(u+v +f]). The transfer matrix t(u, uj) has an interesting 
graphical interpretation ||1 If . Let o: = (ai, . . . , cun ) denote a state in (C") 07V . Then the 
matrix element t(u, ijj)a. is given by 

t(u, u)g = K+(u) (T bc (u, u,))^ K- d {u) (f da (u, , (2.19) 



71 —In 



and is represented graphically by Fig. 6. 
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Figure 6: The Sklyanin transfer matrix t(u,uj). The positioning of the "ra- 
pidities" m ± Ui is crucial. All internal edges (including the ones labelled by 



a, b, c, (i are summed over. 



The commutativity of t(u, u>) implies integrability of the open quantum spin chain 
with Hamiltonian \M 



N-l 



H=J2 H k,k+i + 5 K - (0) + 



k=l 



tr K + (0)H Nfi 
tvK+(0) 



(2.20) 



whose two-site terms are given by 



du 



(2.21) 



u=0 



in the standard fashion.^ However it is not clear at this stage in what sense t(u, ui) 
is a "transfer matrix". One answer is the following | IT| : If the inhomogeneities are 
chosen appropriately then the transfer matrix can be interpreted as that for a vertex 
model on the lattice C More specifically, define the function uj a {u) by 



Kw) ? = h 



i+i 



u, 



(2.22) 



j = 1,...,N. Then t{u 1 u) A {u)) is related to the transfer matrix £d(u) for a vertex 
model on the lattice C. To see this it is most instructive to proceed graphically. Firstly, 
note that regularity of R{u) implies that R(0) has the graphical description given in Fig. 
7. It is then apparent from Fig. 6 that t(u,uj A (u))^ has the graphical interpretation 

3 If the R matrix is in the "homogeneous gauge" so that [Ri2(u), Ri2(v)] = 0, then the 
choice K~(u) = 1, K + (u) = M yields a quantum group-invariant Hamiltonian [^l], p2[ . 
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Figure 7: Graphical depiction of R(0) b J a = p(0) 1/2 5%5 b Q . 
given in Fig. 8 (assuming N is even). Define now the "left and right boundary weights" 

cd 

r b a (u) = K- b (u), (2.23) 
which are depicted graphically in Fig. 9. It is now clear that Fig. 8 can be deformed into 
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Figure 8: Graphical depiction of t(u, u A {u)) 1 the Sklyanin transfer matrix 
with alternating inhomogeneities. iV = 6 in this figure. 

Fig. 2 (modulo an overall scalar factor ^(O)^ 2 ), upon replacing u by u/2. The naming 
of l h a {u) and r b a (u) is also justified even though, strictly speaking, these are boundary 
weights for a model with "spectral parameter" 2u. The left and right boundary weights 
can be separately normalized, as can the bulk weights determined by the i?-matrix. In 
practice, we will find it convenient to choose a normalization such that the bulk weight 
u>i4 and boundary weights and (i.e. where all edges are in the up-arrow state) 
are all unity. Thus we have the key relationship 



*r/o / „ x-JV+i t (u/2,uj A (u/2)) 



(2.24) 
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which defines an integrable vertex model on C The eigenvalues Ad (it) of this transfer 
matrix are obtained once we have diagonalized the Sklyanin transfer matrix t(u, uj). To 
do this, the algebraic or the analytic Bethe ansatz (suitably generalized) will need to 
be used, depending on the model. 

Given an .R-matrix R(u) which satisfies the Yang-Baxter equation, there are typ- 



ically many inequivalent solutions to ( 2.11 ) for K (u) (which are, of course, defined 



only up to an overall scalar factor since ( |2.11|) is a homogeneous equation). In this 



paper we will only be concerned with diagonal .fT-matrices. These lead to boundary 
weights r h a (u) and l h a {u) which are non-zero only if a = b; i.e. to models for which there 
is "conservation of arrows across diagonals" and which can, in principle (but hardly in 
practice) also be solved with the coordinate Bethe ansatz. 




= K a : in) 



Figure 9: Graphical depiction of left and right boundary weights l b a (u) and 
r\{u). 

The "alternating inhomogeneities" uj a {u) also play a similar role in constructing 
diagonal-to-diagonal transfer matrices from inhomogeneous row-to-row ones in the case 
of periodic boundary conditions (cf. Refs. |2^, Q). It is important to appreciate that 
we do not have [tr>(u), tn{v )] = for the chosen w-dependent inhomogeneities. Thus 
we see that "integrability" on the lattice £ is inherited from a more abstract problem 
- commutativity of Sklyanin transfer matrices - instead of being manifestly so. 



3 The six-vertex model 

In Section 2 we described a procedure based on the quantum inverse scattering method 
for obtaining an integrable vertex model on the square lattice with open boundaries, 
related to well-known models integrable on the square lattice with periodic boundaries. 
This procedure has not previously been carried out in its entirety for any vertex model, 
although it has been pushed through [[□]] for the six-vertex model in a special case (with 



one special choice of .fT-matrices) . In fact, solving this model with a general .fT-matrix 
does not involve much more work, especially since the diagonalization of t(u, uj) has 
already been done in Ref. 0. We include it here for completeness and for pedagogical 
reasons. 
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The six-vertex model is associated with the vector (spin-|) representation of the 
(quantum) affine Lie algebra Its i?-matrix reads: 



/ a(u) 



R{u) 



b{u) c{u) 
c{u) b{u) 



\ 



a{u) j 



(3.1) 



where 



a(u) = sinh(w + A), 
b{u) = sinh(-u), 
c{u) = sinh(A), 



(3.2) 



and from this the monodromy matrices T(u, u>) and T(u, ui) are constructed according 
to ( p.2|) and ( p.3|) . This i?-matrix has the properties of regularity and unitarity, with 
p(u) = sinh(A + u) sinh(A — u), P- and T-symmetry separately and crossing-symmetry 
with M = 1 and rj = A. The most general diagonal solution |7| for K~{u) is given 
(up to an unimportant overall scalar function of u) by K~(u) = K(u,£-) where 



sinh(£ + u) 



sinh(£ — u) 



(3.3) 



with £ being a free parameter. Because of the automorphism (|2.13|) , solutions for K + {u) 
are given by K + (u) = K(—u — X, The transfer matrix ( |2.16| ) with these ^-matrices 
has been diagonalized with the algebraic Bethe ansatz, modified for open boundaries (7j . 
A main step in this method involves writing the doubled monodromy matrix U(u, u>) 
in the form 

and using the "intertwining relation" ( [2.15 ) to derive the commutation relations 



B(u)B(v) 
A{u)B{v) 



B(v)B(u), 
sinh(A + v 



u) sinh(w + v) 



B(v)A(u) - 
B(u)A(v) - 



V{u)B{v) 



+ 



+ 



sinh(A +u + v) sinh(v — u 

sinh(u + v) sinh(A) 
sinh(A + u + v) sinh(t> — u] 
sinh(u — v + A) sinh(u +v + 2A) 

sinh(t> — u) sinh(A + u + v) 
sinh(A) sinh(2w + 2A) 



sinh(A) 



sinh(A + u + v 
B{u)V{u) 



-B(u)V(v), 



sinh(t' — u) sinh(2t> + A) 
sinh(A) sinh(2t> ) sinh(2w + 2A) 
sinh(2t> + A) sinh(-u + v + A) 



B(u)V(v) 

B(u)A(v), 



(3.5) 



10 



where T>(u) = sinh(A + 2u)T>(u) — sinh(A)^4(-u). The transfer matrix t(u,uj) can now 
be written (using (p!6|) , (gig]) and Q) as 

sinh(£+ + u + A) ~ sinh(£+ - u) sinh(2w + 2A) 

i u,w = — — — V\u) H — — vAfii). 3.6 

v ' ; sinh(2w + A) v ; sinh(2w + A) y ' v ' 

The algebraic Bethe ansatz says that t(u,uj) is diagonal on the states \u±, . . . ,um) = 
B{ui) . . . B(u M )\Q), where the pseudovacuum \Q) is the state with all arrows up. The 
component B(u) of U(u, uj) acts as a creation operator while C(u), which kills acts 
as an annihilation operator. Furthermore both A(u) and T>(u) are diagonal on \Q). 
To obtain the eigenvalue A(u, uj) of t(u, uj) on \u±, . . . , %), one uses the commutation 
relations (|3.5|) to move A{u) and T>(u) in ( |3.6|) past the B(uj), thus obtaining the 
"wanted terms" 

, , . , . sinh(£ + — u) sinh(2w + 2A) sinhf-u — A — Uj) sinh(w + «,) 
A(u,u>) = a(u,u) . ~~r -[[ 



sinh(2w + A) sinh(w + A + Uj) sinh(w — Uj) 

~ sinh(£ + + u + A) il sinhfu + A — ua) sinhfw + 2A + ua) 

+ o(u,uj) — — — — — - — r — (3.7) 

smh(2u + A) smh(u + A + Uj) smh(w — Uj) 

where a(u,uj) and S(u, uj) are the values taken on by A(u) and V(u) respectively. 
The "unwanted terms" in this procedure cancel (and A (it, uj) becomes analytic in u) 
only if the numbers u k satisfy, for all k, the equations 

a(u k , uj) sinh(u fc + £+) sinh(2w fe ) sinh(zt fc + uj + 2A) sinh(-u fc - uj + A^ ^ 

6(uk, u) sinh(ii fe - + A) J^ k sinh(w fc + Uj) sinh(w fc - Uj - A) 

It is not hard to show - e.g. by graphical means (Fig. 6) - that 

N 

a{u,uj) = sinh(£_ + u) a(u + u>i)a(u — u>i), 

i=i 

N 

6(u,u>) = sinh(£_ — u — A) sinh(2w) \\ b(u + uJi)b{u — uii). (3.9) 

i=l 



After the shift Uj — >• Uj — A/2, the Bethe ansatz equations ( |3.8| ) take on the symmetric 
form 

[ sinh(Mfc — £ + — |) sinh(u fc + — ~) 1 JL sinh(ttfc + u){ + |) sinh(u fc — cjj + 4) 



n 



sinh(w fc + £ + + f ) sinh(u fc - £_ + |) J ^ sinh(w fc + ^ - |) sinh(u fc - ^ - |) 

_ -rj sinh(w A . + Mj + A) sinh(M fc - uj + \\ 
j-J- k sinh(w fc + Uj - A) sinh(« fe - Uj - A) 

The equations ( p.7|) and ( ft.lOQ were obtained in |7j] in the study of the spin-| XXZ 
chain with open boundaries. 
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With the specialization of the inhomogeneities to uj = uj a (u) defined by ( 2. 22] ), the 
second term of A(u,u>) in (|3.7Q vanishes due to 5(u, uj a {u)) = 0. Thus we have (after 



the shift Un 



Uj - A/2) 



A{u, uj a (u)) = sinh(£+ - u) sinh(2w + 2A) sinh(£_ + u)a{2u) N - 1 p(0) N/2 A{u), (3.11) 



where 



. , s ¥L sinh(u — Uj 

a(u) = n 



|) sinh(u + Uj 



(3.12) 



7[ sinh(w + Uj + |) sinh(w — Uj + |) 

From the definitions (|2.23|) we find that the left and right boundary weights are given 
by 



rl{u) 
l\{u) 

ll(u) 



sinh(£_ + u) 

sinh(£_ — u) 

sinh(2w + 2A) sinh(£, 
sinh(£ + + it) 



a 



l\{u) 



sinh(£ H 



u 



(3.13) 



Therefore the eigenvalue A D {u) of the transfer matrix t D {u) is (according to ( |2.24| ) and 

A D (u) = A{u/2), (3.14) 



with A(u) defined by ( ft.12 ). The associated Bethe ansatz equations are given by 



sinh(w fc - £ 4 



sinh(-u fc + £_ 



sinh(w fc + £ + + |) sinh(u fc — £_ + 



sinh(w& + f + |) sinh(-u fc 



u _i_ A N 
2 ' 2/ 



-i N 



sinh(w fc + f — |) sinh(w fc — § — \, } 



M 

n 



sinh(-Ufc + Uj + A) sinh(-Ufc — Uj + A" 
sinh(-u,fc + -Uj — A) sinh(-Ufc — Uj — A, 



■3.15) 



When £± — > — oo, the term in braces on the left hand side of ( |3.15| ) becomes unity.0 
This is the case for which an integrable vertex model on £ has been previously worked 
out [i f] using the same procedure. This case has also been obtained using the coordinate 
Bethe ansatz in Ref. || and is referred to there as the Potts model with free boundaries. 
In fact, the coordinate Bethe ansatz also allows the general case to be solved; the Bethe 
ansatz equations obtained in Ref. H can presumably be recast in the form (3.15p. 



4 With this choice of i^-matrices, the associated open boundary spin-^ XXZ chain 



defined by (2.2C) has J7 9 (su(2))-invariancc El t2U h] 
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4 The Zamolodchikov-Fateev model 



The simplest three-state vertex model is the Zamolodchikov-Fateev 19-vertex |25j or 
the Ai model in the symmetric rank-2 tensor (spin-1) representation [^E] and can 
be constructed from the six-vertex model using the fusion procedure. The associated 
quantum spin chain is the spin-1 XXZ model. The i?-matrix is given byQ 



R{u) 



a 

b 

c 


d 

e 


f 


d 

e 


b 

9 

b 


e 

d 


f 


e 

d 


c 

b 

a 



\ 



(4.1) 



with 



a(u) 
b{u) 
c(u) 
d(u) 
e(u) 

/(«) 
9(u) 



sinh(w + 2A) sinh(w + A), 
sinh(w) sinh(u + A), 
sinh(w) sinh(-u — A), 
sinh(2A) sinh(w + A), 
sinh(w) sinh(2A), 
sinh(A) sinh(2A), 
b(u) + f(u). 



This i?-matrix is regular and unitary, with p{u) = sinh(w + A) sinh(w — A) sinh(w + 
2A) sinh(ii — 2A), P- and T-symmetric and crossing-symmetric with M = 1 and 77 = A. 
The most general diagonal solution for K~(u) has been obtained in Ref. [17] and is 
given by K~(u) = K(u,£-) where 



( sinh(tt + £) sinh(w — A + £) \ 

— sinh(w — £) sinh(w — A + £) 
\ sinh(-u — £) sinh(w + A — £) / 



(4.2) 



Once again, due to the automorphism (|2.13| ) solutions for K + {u) are given by K + {y) = 
K-(-u-\,Z + ). 

The transfer matrix t(u, u>) with these general .fT-matrices but without inhomo- 
geneities Uj has been diagonalized in Ref. JT7J by generalizing the approach used to 

5 This is also the i?-matrix associated with the vector representation of 
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solve the corresponding model with periodic boundaries. The model with inhomo- 
geneities can be treated likewise. This involves the construction of an auxiliary transfer 
matrix 

J(u, w) = tr (^K + {u) c T{u) a K-{u)^f{-u)) . (4.3) 

which commutes with t(u, u>). The auxiliary transfer matrix corresponds to that for a 
model whose monodromy matrices are constructed from a simpler .R-matrix a R(u) and 
whose i^-matrices a K ± (u) are those from the six- vertex model (|3.3|) , but with the free 



parameters shifted according to 
in turn is given by 



- | and £ H 



-£ + + w. The i?-matrix a R(u) 



a 








b 




d 






C 




d 


d 




c 






d 




b 








a 



(4.4) 



with 



a 
b 
c 
d 



sinh(w + 4p), 
sinh(-u + |), 
sinh(w — !), 
[sinh(A)sinh(2A)] 1/2 



and is associated to the tensor product of a spin-i and a spin-1 representation of A\ 



(i) 



The monodromy matrices a T(u) are 2x2 and thus the auxiliary transfer matrix 
a t(u, uj) can be diagonalized with the algebraic Bethe ansatz |]T7[ in the same manner 
as for the six- vertex model in Section 3. Leaving out the details, we find that the 
eigenvalue in the inhomogeneous case is given by 



where 



r A(u,w) 



smbj2^ + 2A) N 
smh(2w + A) 



Ui)A(u) 



+ P(u) 



sinh(2tt) 
sinh(2w + A) 



N 



Yl c(u + uJi)c(u - uJi)B(u) 



a{u) 


= sin 


P{u) 


= sin 




M 


A{u) 


= n 




3=1 




M 


B(u) 


= n 




3=1 



2 

3A 



+ m sinh(w + + 4) 



sinh(w 



Mi 



sinh(u + Uj 



\ sinh(u — Uj + |) sinh(-u + Wj + | 



sinh(-u 



Ui 



+ ^)sinh(w + M i + 4r) 



sinh(w — + |) sinh(-u + Wj + |) 



(4.5) 

(4.6) 
(4.7) 

(4.8) 
(4.9) 
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The eigenvalue for the original transfer matrix t(u, u>) is related to that for the 
auxiliary one by 

A(w, u>) = 7(u) ff A(w - |, u>) a A(u + |, w) + <J(u) (4.10) 

where we have defined 

sinh(2w) sinh(2w + 2A) 



7(u) 



sinh(2w + A) 



2 



. . , . sinh(2w + 3A) sinh(2w - A) SL _. ._, . 

* « = -a(«- 1/50+1 ■ r J - , 2 II < u + <*) e (« ~ 

smh(2w + X) z " 

e(u) = sinh(w + 2 A) sinh(w - A). (4.11) 

The Bethe ansatz equations satisfied by the roots Uj are the same ones as for the 
auxiliary transfer matrix, namely for all k: 

sinh(M fc + - A) sinh(w fc - £+) £L sinh(u fc - u { + X) sinh(ii fc + uii + A) 
sinh(u fc - + A) sinh(u fc + £+) ^ sinh(tt fc - - A) sinh(u£; + w» — A) 

-rj sinhpfc - ttj + A) sinh(w fc + uj + A) 
^ sinh(w fe - - A) sinh(w fc + uj - A) ' 

When co>j = these are the Bethe ansatz equations which determine the energy spectrum 
for the spin-1 XXZ chain with open boundaries \T%. If, in addition, we choose £± — > — oo 
then the spin chain has Lf 9 (sw(2))-invariance [|i~7 |. 

With the alternating inhomogeneities lj = lj a (u), the transfer matrix eigenvalue 
evaluates to 

A(u,u A (u)) = a(w-|)a(w + |)sinh(2w + 2A)sinh(2w + 3A) 

x a(2u) N - 1 p(0) N ^A(u). (4.13) 

Once again, only one term survives. This seems to be a generic feature of all ver- 
tex models on the lattice C, as the coordinate Bethe ansatz would also confirm. By 
substituting the expressions for the R- and i^-matrices into ([2.231) , the left and right 
boundaries l b a {u) and r b a {u) are found to be 

l\(u) = sinh(2u + 2A) sinh(2u + 3A) sinh(£ + — u) sinh(£ + — A — u), 
,2/ \ ii ( N sinh(^ + - X + u) 



ll(u) = l\{u 



sinh(£ + — A — u) ' 
sinh(£ + + u) sinh(£ + — X + u) 
sinh(£ + — u) sinh(£ + — A — u) 
sinh(u + £_) sinh(u — A + £_), 



2 x sinh(£_ -u) 

r 2\ u ) = r ii u ) 



sinh(-u + £-) ' 
r 3 (u) - r i fu) sinh (^-e-)smh(» + A-e-) 
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It is then straightforward to show that l\(u)r\(u) = sinh(2it + 2A) sinh(2w + 3X)a(u — 
T?)a(u+ij). Thus the eigenvalue expression for the diagonal-to-diagonal transfer matrix 



to{u) is simply (cf. equations Q4.13|) and (|2.24|) ) 

A D (u) = A(u/2), 



(4.15) 



where A(u) is defined in 



if we normalize all Boltzmann weights such that W14 



■w 



,R 



w{" = 1. The corresponding Bethe ansatz equations are given by 



sinh(w fc + — A) sinh(-u fc — 
sinh(M fc - £_ + A) sinh(w fc + 



M 

n 



sinh(w fc — | + A) sinh(w fc + | + A) 
sinh(-Ufc — | — A) sinh(-u fc + | — A) 

sinh(-Ujt — Mj + A) sinh(wfc + ttj + A) 



(4.16) 



The case where £_ 
the coordinate Bethe ansatz |28l . 



sinh(u fc - Uj - A) sinh(« fe + Uj - A) 
u/2 + 2A and £ + = —u/2 — A has previously been obtained via 



5 The A { 2 ] vertex model 

The i?-matrix associated with the affine algebra A% in the vector representation is 
given by |2T 



where 



R{u) 



a 

b 

b 


c 


> 

c 


d 


b 

a 

b 


c 


d 


d 


b 

b 

a 



\ 



a[u) 
b(u) 
c(u) 
d{u) 



sinh(-u + A), 
sinh(-u), 
e u sinh(A), 
e~ u sinh(A). 



(5.1) 



This i?-matrix has the properties of regularity and unitarity, with p(u) = sinh(A + 
u) sinh(A — u), PT-symmetry and a weak version of crossing-symmetry |TJ], ^TJ 

-i - o 



(Rl2(u) 



t2 



f(u) M Ri 2 (u + 2rj) M~\ 



(5.2) 
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with 77 = 3 A/2, M 



diag(e 2A , 1, e 2A ) and scalar factor 
sinh(tt) sinh(u + 3A) 



f(u) 



sinh(M + 2A) sinh(u + 4A) 



(5.3) 



With the property (|5.2|) replacing crossing-symmetry, the condition fl2.12|) on K + {u) 
for integrability remains valid |H| , with the automorphism ( 2.13| ) again holding. It is 
easy to see that K~(u) = 1 is a solution of (|2.11|) for the present model, as with all the 
models considered here. The most general diagonal solution |RTJ for K~(u) is given by 
K~{u) = K(u,£,_,l_) where 



sinh(£ — u)e u 
sinh(£ + u)e~ u 



\<a<l 
I + 1 < a < 3 



(5.4) 



where £ is arbitrary and I can be either 1, 2 or 3. By ( [2.13] ) solutions for K + (u) follow, 
with corresponding £ + and l + . 

The most general commuting transfer matrix t(u,uj) then depends on the four pa- 
rameters £±,^±- This transfer matrix has recently been diagonalized with the nested 
Bethe ansatz generalizing an earlier result Jl8| valid for the case K~{u) = 1 and 
K + (u) = M, which leads to an [7 g (sM(3))-invariant spin chain ||21|| . In fact these authors 
studied the whole family of A^\n > 1 vertex models. We will only be interested here 
in the three cases where /_ = /+. Note that the case /+ = /_ = 3 is actually equivalent 
to the [/g(sw(3))-invariant case, by rescaling the i^-matrices. The expression for the 
eigenvalue A(u, u>) can be written as 



N 



A(u,u) 



a(u) 



where 



+ P(u) 
+ l(u) 

A(u) 
B{u) 
C{u) 
D(u) 



sinh(2w + 3A) 
sinh(2w + A) ■ , 

sinh(2u) sinh(2u + 3A) 



a{u + u>i)a(u — Ui)A(u) 

N 

Yl b(u + Uj)b(u - Uj)B(u)C(u) 



sinh(2w + A) sinh(2w + 2A) 
sinh(2u' 



N 



sinh(2u + 2A) 



Y[ b(u + u)i)b{u — uJi)D{u) 



M 



„ sinh(w + fif' 
II 



sinh(w 



j=i sinh(-u + fij + |) sinh(w 



M 

n 



2 

j. sinh(w + fj,) l) + 4r) sinh(w - + ^ 



1 sinh(ii + /iy + |) sinh(w — /j,f ' + |) 



M 2 



sinh(w + + |) sinh(w 



n 

i=i 

n 



u (2) 



(2) 



j=i sinh(w 
M2 sinh(w + 



+ f ) sinh(w - nf + |) 

5A 



+ tt) sinh(w — ji 



(2) 



3A^ 
2 / 



1 sinh(u + jjf 3 + 4p) sinh(ii — (if' + 



3V 



4 



(5.5) 

(5.6) 
(5.7) 
(5.8) 
(5.9) 
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The (nested) Bethe ansatz equations follow from analyticity of A(u, oS) and are given 
by 



-|)tt sinh(/if + Wi + |) sinh(/if - ^ + §) 



/^f - |) *= A i sinh(^ i; + w< - |) sinh(/i^ i; - ^ 



n 



sinh(/i^ ; + /if + A) sinh(/i^ ; - /if + A) 



^ sinh(/if + /if — A) sinh(/if — /4 — A) 
x sinM/xf + /if - A) smru>g - /if - |) 
M sinh^f + /if + |) sinh(/if - /if + |: 
- |) S sinh(/if + /if + A) sinh(/if - /if + A) 



n 



l(^k' - I) &k sinh(/if + /if - A) sinh(/if - /if - A) 

x jr sinh(/ij 2) + /if - |) sinh(/if - /if - |) 
X M sinh(/if + /if + |) sinh(/zf - /if + |) 



for all k. For the three types of boundaries under consideration the scalar factors a, (3 
and 7 are given by 
m l_ = U = 1: 



= e2 sinh(£_ — -u) sinh(£ + + u 



A 



/?(u)=7(ii) = e2 sinh(£„ + w + A) sinh(£+ - u - (5.12) 



(ii) L = l + = 2: 



A 



Qf(«) = P(u) = e 2 sinh(£_ — u) sinh(£ + + u + 

A 



7 ( M ) = e ~ 2 s inh(£_ + it + 2A) sinh(£+ - w - §p) (5.13) 
(iii) /_ = /+ = 3: 

_3A 

a(w) = /3(«) = 7(w) = e 2 sinh(£_ - w) sinh(£ + + u + ^). (5-14) 

These expressions for a(w), /3(w) and 7(14) follow from the results in Ref. [|1^], after 
a shift in £ + necessary to match up with the different parametrization of K + {u) used 
therein. 

We now specialize the inhomogeneities oj to u> (u). As usual, only the first term in 
the eigenvalue expression (|5.5| ) survives: 

A(u, u A {u)) = a{u) sinh(2w + 3X)a(2u) N - 1 p(0) N/2 A(u). (5.15) 

The right and left boundary weights r b a {u) and l h Ju) follow from the definitions (|2.23|) , 
and for the three cases of interest read: 
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(i) L = L 



1: 



(ii) L 



r 2 2 (u) ■■ 



l\{u) 

= ll( U ) 

r\{u) 
rf(u) 



sinh(2w + 3A) sinh(-u + 

Me 2 ", 

sinh(£_ - tt)e", 



r\[v. ; 



( sinh(£_ + it) 
sinh(£_ — u) 



-2u 



(5.16) 



l\{u) = ll{u) 



r\{u) 



rf(u) 
rf(u) 



sinh(2u + 3A) sinh(u + £+ + |)e- A/2 -", 



sinh(§ + £+ + u) 
sinh(£_ - u)e u , 
1 sinh(£_ + u) 2u 



sinh(£_ — u) 



(5.17) 



fiii) I- = L 



3: 



l\{u) = lj(u) = ll(u) 

r l( u ) = r l( u ) = r i( u ) 



sinh(2u + 3A) sinh(u + + ^)e~ 
sinh(£_ -u)e u . 



(5.18) 



From these we find that l\(u)r\(u) = a(u) sinh(2-u + 3 A) for all cases. Therefore if we 
normalize all the vertex weights such that wu = 1 = wf = lyf, the eigenvalue for 
the diagonal-to-diagonal transfer matrix to{u) becomes (cf. (|2.24j) and ( |5.15|) ) simply 
Ad(m) = A with A(u) defined in fl5.6|) . The corresponding Bethe ansatz equations 
are given by (|5.11 ) and 



"Of 



(i) 



sinh(/i^ 1) + f + |) smh(^ ] 



- + 

2 ~ 2 



|) [sinh(/if + |-|) sinh(/if - | - |; 
Ml -"M/if + /if + A) sinh(/if - /if 



N 



n 



y-* sinh(/i^ 1) + /if - A) sinh(^^ - /if - A) 



x 



Ma sinh(/if + //f - |) sinh^f - (2) 



n 



.y— 1 smh(4 1} 



M 2) 



) sinh(/if - /if + |) 



(2) 



(5.19) 
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6 The A? vertex model 



The i?-matrix for the Izergin-Korepin p2[ or A^p model [IR [L6| is given by 



where 





( ° 


b 

d 


e 

9 




f 






e 


b 






R(u) = 




9 


a 


b 


9 

e 






I 


9 


e 


d 

b 












c 



\ 



(6.1) 



a(u) 
b{u) 
c{u) 
d{u) 
e(u) 
e(u) 

m 

/(«) 
9(u) 
9{u) 



sinh(u — 3A) — sinh(5A) + sinh(3A) + sinh(A), 
sinh(u — 3A) + sinh(3A), 
sinh(-u — 5A) + sinh(A), 
sinh(u — A) + sinh(A), 
-2e~5 s inh(2A) cosh(f - 3A), 

u 

-2e2 sinh(2A) cosh(| - 3A), 

-2e- u+2X sinh(A) sinh(2A) - e~ x sinh(4A), 

A 



„u-2A , 
-^+2A, 



2e u ~ ZA sinh(A) sinh(2A) - e A sinh(4A), 



2e"2 +ZA sinh(|)sinh(2A), 
-2e2" 2A sinh(|)sinh(2A). 



This .R-matrix has the properties of regularity and unitarity, with p(u) = [sinh(A) + 
sinh(w — 5A)][sinh(A) — sinh(u + 5A)]. It is also PT-symmetric and crossing- symmetric, 
with 7] = — 6A — m and 

M = diag(e 2A ,l,e~ 2A ). (6.2) 

Diagonal solutions for K~{u) have been obtained in |T3|. It turns out that there are 
only three solutions, being K~{u) = 1, K~{u) = T + (u) and K~{u) = T~(u), with 



[u 



0» 



(6.3) 
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where we have defined 



^(u) = cosh(§ -3A) ±isinh(§) (6.4) 
0±( M ) = cosh(^ + 3A) ^isinh(^). (6.5) 

There is an absence of free parameters in these solutions, in contrast to the previous 
models. By the automorphism (|2.13|) , three solutions for K + (u) follow. We have thus 
nine possibilities for the commuting transfer matrix t(u, uj). We will only be interested 
in the three cases 

(i) K~{u) = 1, K+(u) = M, 

(ii) K~(u) = T+(u), K + {u) = T + (-u + 6A + mfM, 

(iii) K-{u) = V-(u), K+{u) = T-(-u + 6A + mfM. 

The transfer matrix for case (i), whose corresponding spin-chain is {7 g (sM(2))-invariant 



21], has been diagonalized by the analytic Bethe ansatz in Ref. |33[ for the homoge- 



neous case LUi — 0. It is not difficult to introduce the inhomogeneities, and the result 
for the eigenvalue can be written as 

a/ \ tt / \ / \ sinh(w — 6A) coshfw — A) ., . 

a(«,«) = n + ->m* - "0 sinh( „ _ 2A / C03h( ; _ 3A ) mu) 

smn(« — 2A) smn(« — 4A) 

tt ,/ n ,/ \ sinh(w) cosh(w — 5A) . . 

+ Y\ d(u + uAd(u - Ui) — — ~r~\~\ \ 7 — — -C7(it) , 6.6 

^ v i; v * ; sinh(w-4A)cosh(w-3A) K h K ' 

with 

'' sinh(±(n + m 3 -) + A) sinh(|(w - %) + A) 
W | = i 1 sinh(i( M + Mj ) - A)sinh(i(u- Uj ) - A) 1 J 

. -pi sinh(i(w + ttj) — 3A) sinh(i(u — Uj) — 3A) 

j-J^ sinh(i(u + Uj) — A) sinh(i(u — Uj) — A) 

x y\ cosh(±(u + Uj)) cosh(±(u - Uj)) 

^ cosh(±(w + Uj-) - 2A) cosh(±(w - Uj ) - 2A) 1 ' ' 

"' cosh(±( M + ^) - 4A) cosh(±( M - Uj) - 4A) 
/=i cosh(i(w + Ui ) - 2A) cosh(i(w - Uj ) - 2A) ' 1 ' ' 

For the above case, quantum group invariance was used to determine properties 
of the transfer matrix eigenvalue essential for the analytic Bethe ansatz to be pushed 
through. For the non-quantum group invariant cases, the analogous calculations have 
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not been carried out. Nevertheless, we expect the eigenvalue expression to be simply 
modified to 

N N 

A(it, u>) = FT c(u + uji)c(u — Ui)at(u)A(u) + TT b(u + 0Ji)b{u — uJi)(3{v)B{v) 

i=l i=l 
N 

+ Hdiu + Uijdiu-uJi^iujCiu), (6.10) 
i=i 



with A(u), B(u) and C{u) defined as in (|6.7|) , (|6.8|) and ( |6.9|) and the boundary- 
dependent terms a(u), (3{u) and 7(11) are to be determined by the action of t(u,u>) 
on the pseudovacuum \Q) (state with all arrows up). By a straightforward calculation 
- e.g. by graphical means (Fig. 6) - we find them to be given by 

„ K+ 2 e 2 K+ f- 2 g 2 e 2 \\ 



c 2 — b 2 c 2 — d 2 V c 2 — b 2 I \ 
K+K^e 2 K+K^g 2 K+K u g 2 e 2 



13 K ™ K ™ C 2_ b 2 d 2_ b 2 ( C 2_ 6 2)( 6 2_ rf 2) 

7 = K ^ K ^ + d^-^dA f -w^i 2 "- (6 - n) 



Naturally, for case (i) the result (|6.6j ) is recovered. For the other cases ( |6.11|) yields 
(ii): 

(ip + ) 2 

a ( u ) = ^7TT^ + ( M ) sinn ( M_6A )' 
c(2u) 

sinh(w) cosh(A) + i sinh(2A — u) 

p(u) = a(u) 



sinh(w — 4A) cosh(A) — isinh(2A — u) 



in : 



g ju) sinh(w) sinh(M - 2A) 

7(w) = «(«)——— — — — — — , (0.12) 

,K ' v '£ + {u) sinh(u - 6A) sinh(u - 4A) ' v ; 



(VO 2 

= tttt^ (m) smh(w - 6A), 

sinh(-u) cosh(A) — i sinh(2A — u) 

p(u) = a[u) 



sinh(w — 4A) cosh(A) + i sinh(2A — u) ' 

£+(n) sinh(w) sinh(n - 2A) 

7(m) = a(u) - , . — :— n ttt, (o.ld) 

' V ; v ^-(w)sinh(w-6A)sinh(w-4A)' v ; 

where we have defined 

^(u) = 2 (cosh(3A -u)±i sinh(2A)) . (6.14) 
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The Bethe ansatz equations associated with the eigenvalue expression (|6.10|) are given 
by 



n 

2=1 



c(u k + 2A + uJi)c(u k + 2X- lu 1 
b(uk + 2A + uji)b{uk + 2\ — u)i 

smh[\{u k 



a{uk + 2A) sinh(ufc + 2A) cosh(tifc — A) 
P(uk + 2A) sinh(ufc — 2A) cosh(wfc + A) 



Uj 



2 A] sinhf^fM/c 



2A] 



sinh[i(w fc + uj) + 2A] sinh[|(u fc - Uj) + 2A] 
x cosh[i(n fc + Uj) + A] cosh^Qfc - Uj) + A] 
cosh[|(ttfe + Uj) — A] cosh[!(wfc — Uj) — A] 

Following the general procedure laid down in Section 2, we set uj = uj a (u) in flOOP . 
As usual, only the first term survives: 

A(u, u; A (u)) = c(2w) Jv c(0) Ar Q;(w)^(w). (6.16) 

From the definitions (p.23|) the left and right boundary weights l h a (u) and r b a {u) are 

found to be given by 

(i): 



l\{u)= ll{u) = l\{u) 
r\(u) =rl(u) =rl(u) 



2 sinh(-u — 6A) cosh(w — A) 
1, 



(6.17) 



(iii): 



l 2 2 (u) 

l 3 3 (u) 
r\{u) 

r 2 2 {u) 

r\{u) 

l\{u) 
lliu) 
lliu) 



r}(u) 



r\(u) 
rf(«) 



i) + (u) sinh(M-6A)eVH 

ip + (u) 
l\(u)e- 2u 

ij + (u)e~ u 

%p + {u) 
r{(u)e 2u , 



-ij) (u) sinh(-u — 6A)e M £ (u) 



u 



(6.18) 



ip (u) 
l\(u)e- 2u 

ip-(u)e~ u 



[u 



r l(u)- , 
V ( ) 

r\(u)e 2u . 



(6.19) 
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For all three cases, we have l\(u)rl(u) = c(2u)a(u). Hence, by normalizing all bulk 
weights such that Wu = 1 and boundary weights such that u>f = u>f = 1, the eigenvalue 
for the diagonal-to-diagonal transfer matrix tr>(u) is simply 

A D (u) = A (|) , (6.20) 

where A(u) is defined in (|6.7|) . The corresponding Bethe ansatz equations are given by 
(|6.15|) with the left hand side replaced by 



c(w fe + 2A + f)c(w fc + 2A-§) 



b(u k + 2A + f)b{u k + 2A - f 



2/J 



a(u k + 2A) sinh(w fc + 2A) cosh(w fc — A) 
f3{u k + 2A) sinh(wfc — 2A) cosh(ufc + A) 



(6.21) 



For case (i), the term in braces in (|6.21| ) becomes unity (cf. ( |6.6| )), whereas for cases 
(ii) and (iii) it simplifies to 



/ cosh[i(-» fc - A + f)] ' 
Vcosh[i(n fc + A-f)] 

and 



(6.22) 



/ sinh[lK-A + f)] V 

Umh[iK+A-f)]; [b - Z6) 



respectively. 



7 Loop models with integrable boundaries 

Let C be the lattice depicted in Fig. 1. A loop model on £ has partition function 

Z Xoop = Y,PT---pTrn p (7.1) 

G 

where the sum is over all configurations Q of non-intersecting closed loops which cover 
some (or none) of the edges of C. The possible configurations at each vertex are shown 
in Fig. 10, with the one of type i carrying a Boltzmann weight p^. In the configuration 
Q, rrii is the number of occurences of the vertex of type i while P is the total number 
of closed loops of fugacity n. Much work has been done on models with periodic 
boundaries (see e.g. Refs. |2y, |4j]), where the boundary weights p w to pi 3 are absent. 
In particular an 0{n) vector model can be mapped onto such a loop model, with n 
being precisely the loop fugacity in Z loop . The n — > limit of such an 0{n) model is 



especially interesting because of its polymer interpretation |35 . 

Instances where the loop model admits exact solutions are known. Such important 
cases can be related to more well-known integrable vertex models via the loop model- 



to- vertex model mapping pi EM. This involves first assigning orientations to the loops 
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Figure 10: Allowed vertices for the loop model with partition function Z\ oop . 
Vertex % carries Boltzmann weight p%. 

in Q . The explicit closed loop counting in Z\ oopi which is a non-local procedure, can be 
eliminated by introducing a weight s for each left (right) turn along an oriented 

loop if s is chosen so that n = s 4 + s -4 . The loop model partition function then becomes 
equivalent to a (three-state) vertex model partition function (|1.1|) , whose weights are 
functions of pi, s and, in general, gauge factors a, b, c and d (whose presence leave the 
partition function unchanged). The relations between the vertex weights Wi, wf and 
wf- and the loop weights p, are given explicitly by (refer to Figs. 3, 4 and 10) 

w\ = Pi Wi = Pw^jh 



(7.2) 



(oc) 

w 2 = p2s{ac) W2 = pn 

1 

^3 = p%s (bd) W3 = p 10 s( bd ) 

: : wf = Pi2s(ac) 

Wl8 = (P8S 2 + P 9 S~ 2 ) f d = P13 

wis = (Pss~ 2 + P9S 2 ) g w§ = p 12 s-\bd). 

The question we want to pose is the following: Given an integrable vertex model on C, 
can a loop model be found? For this to be possible, the vertex weights must satisfy a set 
of constraints implied by the equations ( |7.2j ). The constraints on the bulk weights are 



precisely those for the model with periodic boundaries [^(], [34]]. For open boundaries, 
we see that the constraints 



W3 : w[ = Wi : u>3 = w 2 : w 3 (7.3) 

on the boundary weights must also be satisfied. The important thing to note about the 
criterion (|7.3|) is that it is gauge-independent. 
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7.1 The 0(n) loop model 

To see how this works in practice, we look at the 0(n) model in some detail. This loop 
model is related to the A 2 model |20| . Following Ref. [[34]] we use the parametrization 
of the A 2 .R-matrix in terms of u and A, which is related to the parametrization used 
in Section 6 by u = 2iu and A = iA + m/2. With a suitable choice of gauge factors 
a, b, c, d the bulk loop model weights are given by 





Pi 


= sin(3A — u) sin(£t) + sin(2A) sin (3 A), 


P2 


= P3 


= ei sin(3A — u) sin(2A), 


Pa 


= P5 


= e2 sin(w) sin(2A), 


P6 


= P7 


= sin(3A — u) sin(tt), 




Ps 


= sin(3A — u) sin(2A — u), 




P9 


= — sin(A — u) sin(u), 



(7.4) 

where ef = e| = 1 and n = — 2cos(4A). Referring back to equations ( |6.17| ) to ( |6.19| ) we 
see that for both choices of non quantum algebra invariant boundaries (cases (ii) and 
(iii), but not case (i)) the constraints ( |7.3j ) on the boundary weights are satisfied, the 
ratio being 1 : e u in both cases. Thus, for both sets of integrable boundaries a loop 
model interpretation is available. For case (ii), with the boundary vertex weights given 
by (|6.18| ), the corresponding boundary loop weights are given by 

Pw = pi2 = sin[i(3A -u)], 

Pu = P13 = ei sin[i(3A + u)}, (7.5) 

with the gauge factors used being the same ones as those which led to ( |7.4| ). As with 
vertex models, the bulk weights (pi, . . . , pg), left boundary weights (pio,Pn) and right 
boundary weights (pi 2 , P13) can be separately normalized without changing the essential 
properties of the model. We have chosen one suitable normalization here. 

The relevant Bethe ansatz solution can be lifted from Section 6 and, up to normal- 
ization factors, reads 

A D (u) = f] sinh [l(^ + u i) + jj)] gmh[|(ig - u 3 ) + iA)] ^ ^ 

, =1 sinh[i(in + Uj) — iA)] sinh[i(i-£t — Uj) — iA)] 

where the numbers Uj follow as roots of the equations (1 < k < M) 

( sinh[^(nfc — iA)] \ / sinh[^(Mfe — 2iA — iu)) sinh[^(ufc — 2iA + iu)] \ 
\sinh.[^(uk + iA)] / \sinh[i(-u fc + 2iA + iu)] sinh[i(-u fc + 2iA — iu)] J 

sinh[i(w fc + uj) — 2iA] sinhf^Mjt — Uj) — 2iA] 
sinh[^(Mfc + Uj) + 2iA] sinh[^(wfc — Uj) + 2iA] 

sinh[i(wfc + Uj) + iA] sinh[i(uA; — Uj) + iA] 
sinh[|(-u fc + Uj) — iA] sinh[i(-u fc — Uj) — iA] 



(7.7) 
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A special case of this solution has been derived using the coordinate Bethe ansatz |J. 
This pertains to the case u — A which is equivalent to the 0(n) model on the honeycomb 
lattice ]36], |20], |37]], on account of the vanishing of pg which enables all the other loop 



vertices to be "pulled apart" horizontally. The vertex weights given in Ref. || differ 
from the ones obtained in Section 6 by a gauge transformation. 
For case (iii), the boundary weights ( |7.5|) are replaced by 



Pio = P12 = cos[|(3A-«)], 

Pn = P13 = eicos[i(3A + tt)]. (7.8) 



The corresponding Bethe ansatz solution is given by ( [7.6| ) and (|7.7| ), with the sinh 
function in the squared term on the left hand side of the latter replaced by cosh. The 
honeycomb lattice Bethe ansatz solution corresponding to case (ii) with u = A has been 
analysed for its critical behaviour in Ref. |J in connection with self-avoiding walks. 
A study of the similarlyly interesting honeycomb lattice solution corresponding to case 
(iii) will be presented elsewhere, together with that for the square lattice solutions 
(general u). 

7.2 The A { 2 } loop model 

We now turn to the loop model connected with the A% vertex model |37], |34j|. With a 
suitable choice of gauge factors, the bulk loop weights are given by 

p!=p 8 = sinh(w + A), 

P2=P3 = esinh(A), 

p 4 = p5 = 0, 

Pe = p7 = p9 = sinh(w), (7.9) 

where e 2 = 1 and the loop fugacity is n = — 2cosh(A). As seen from the i?-matrix, 
the ratio W2 : for this model is 1 : e~ 2u . An inspection of equations (|5.16 ) to ( 5.18|) 



shows that the criterion ( |7.3| ) is satisfied only for case (ii) and only if the parameters 
£± are taken to — oo. This represents the only choice of integrable boundary weights 
found in Section 5 for which the model has a loop model interpretation. The boundary 
loop weights work out (using equations ( |7.2| )) to be 



Pio = Pi3 = e u/2 , 

Pn=Pi2 = ee""/ 2 . (7.10) 
The corresponding Bethe ansatz solution reads: 

A D (u) = n ( 7 .n) 

/isinhd + ^+Dsinhd-^ + l) 
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with the Bethe ansatz roots ^ and satisfying the nested equations 



(2) 



sinh(/x[, 1) + | + |) sinh(/4; 



(i) 



u I A" 
2 "r 2> 



sinh^ + 



I) sinh(/i« 



N 



Ml Binh^ + //f ) + A) sinh(//i 1) 



n 



A) 



sinh^ + — A) sinh(/4 



A) 



X 



Ma sinh(//J 1) 1 - (2) 



n 



sinh(/i^ 1) - //V 



,(2) 



i sinh(/4^ + /i^ 

(2) 



sinh(/i^ 



4 2) + |) 



Ma sinh(/il 2) + nf + A) sinh(/i^ 



n 



(2) 



(2) 

Mi 



A) 



"fe sinh(/i^ 2) + /if } - A) sinh(/4 2) - - A) 



x 



g sinh(/.i 2) + ^ 



|) sinh(/4* 



(2) 



2 / 



i= i sinh(/i^ 2) + nf> + |) sinh(/^ fc 



(2) 



^ + I) 



(7.12) 



The honeycomb limit u = —A of the A 2 model has been shown ||38| to yield the 
fully-packed loop model J39| on the honeycomb lattice, which is equivalent to the zero- 
temperature limit of the 0(n) model. For periodic boundary conditions the relevant 
Bethe ansatz solution has been analysed |38| and applied to the problem of Hamiltonian 
walks (i.e. self-avoiding walks which visit every vertex) on the honeycomb lattice. The 
Bethe ansatz solution obtained here for open boundary conditions is relevant to the 
surface critical behaviour of Hamiltonian walks and will be discussed elsewhere. 



7.3 The Zamolodchikov-Fateev loop model 

We will now investigate the loop model connected with the Zamolodchikov-Fateev ver- 
tex model of Section 4. It turns out that a loop model interpretation is possible if and 
only if the loops are allowed to intersect p(| . To take this into account an extra weight 
p (see Fig. 11) is included in the definition of the partition function Z\ oov . Further- 
more, the loop fugacity is fixed at n = 2 and the resulting model has an interpretation 
as the high-temperature expansion of an 0(2) model on the square lattice The 
bulk loop model weights are then found to be 





Pi 


= sinh(w) sinh(u + A) + sinh(A) sinh(2A), 


Pi 


= P3 


= ei sinh(2A) sinh(w + A), 


Pi 


= P5 


= e 2 sinh(w) sinh(2A), 


Pe 


= P7 


= sinh(u) sinh(w + A), 




Ps 


= cosh(-u) sinh(2A) sinh(w + A), 




P9 


= — sinh(w) sinh(2u) cosh(u + A), 




Po 


= sinh(u) cosh(2A) sinh(u + A), 
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Figure 11: Extra vertex with Boltzmann weight po needed if loops are allowed 
to intersect 



where 



1. From the equations ( f4.14p we see that the criterion ( |7.3| ) is satisfied 



and thus a loop model interpretation is possible if and only if we set £ + = = A/2 in 
the integrable weights found in Section 4.[] The boundary loop model weights are then 
found to be 



Pio = Pi2 = sinh[i(w + A)], 
Pn = P13 = eisinh[|(A - u)}. 



The corresponding Bethe ansatz solution is given by 
A D (u) 



jL sinh(| — Uj — I) sinh(| + Uj — ^ 
sinh(| - Uj + §) sinh(| + Uj + I 



with the Uj being roots of the equations 



/ sinh(n fc - |) 
\sinhOfc + |) 



sinh(-u fc — | + A) sinh(w fc + | + A) 
sinh(-u fc — | — A) sinh(-Ufc + | — A) 

sinh(u fc — Uj + A) sinh(u fc + Uj + A) 
Uj — A) sinhtwfc + Uj — A) 



M 



Jj- k sinhOfc 



(7.14) 



(7.15) 



(7.16) 



7.4 Loop models related to the six-vertex model 

To conclude the study of loop models on the lattice C, we now return to the six-vertex 
model. The simplest loop model connected with the six-vertex model is the so-called 
dense (i.e. all edges of £ are covered) 0(n) model with partition function 



J dense 



E 

Q 



mg m 9 mio m l2 P 
J 8 H9 PlO Pl2 a i 



(7.17) 



6 Note that as with the A% and A% models, the choice of if-matrices required for a 
loop model interpretation is not the choice which gives rise to quantum algebra-invariant 
spin chains. 
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and Boltzmann weights 



p 8 = sinh(ii + A) 

p 9 = sinh(tt). (7-18) 

The loop fugacity is n = — 2cosh(A). Since the vertex weights w 2 and w 3 are zero, 
we need to modify the criterion ( [7.3| ) which the boundary vertex model weights need 
to satisfy in order to have a loop model interpretation. It is easily seen that the new 
requirement is given by 

: wf = u>f : = W14 + s 4 u>i 6 : Wig. (7-19) 

The required ratio is 1 : e eu where e 2 = 1, and is achieved only if the free parameters in 
the integrable six- vertex model boundary weights found in Section 3 take on the values 
£± — >• oo or £± — ► — oo.[] The boundary weights can be conveniently normalized to 

pio=Pi2 = l. (7.20) 

The related Bethe ansatz solution is given by 

A ( u \ = TT sinhd-^-^sinhd + ^-A) 
D{ ' f} x sinh(f + Ui + §) sinh(f - u { + |) 1 ' ; 

with Mj being roots of the equations 

smh(uk + § + |) sinh(-u fc — | + |) ^ sinh(-u fc + itj + A) sinh(w fc — Uj + A) 

sinh(-u fc + f — |) sinh(-u fc — | — |) r^, sinh(w fc + Uj — A) sinh(w fc — Uj — A) 

(7.22) 

There is another loop model connected with the six-vertex model, a dilute (i.e. not 
all edges of C are covered) 0{n) model. This model has a partition function Zdiiute 
equivalent to -Zdense but involves different loop weights. The equivalence can be shown 
[UJ] by mapping Z dense to a 2-colour dense loop model |34[] with fugacities ri\ and n 2 
such that rii + ri2 = n. Setting n 2 = 1 and then summing over the second colour one 
arrives at ^dilute with weights 

Pi = sinh(ii + A) + sinh(w) 
P2= p3 = Ps = sinh(u + A) 
Pa = P5 = P9 = sinh(u) 
Pio = Pn = P12 = pia = !• (7-23) 



7 Recall that the choice £± —> — oo leads to a C/ 9 (su(2))-invariant spin-i XXZ chain. 
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The Bethe ansatz solution is the same as for the dense model, but the loop fugacity n 
is now given by n = — 2cosh(A) + 1. This solution has recently been studied in 



in 



relation to walks on the Manhattan lattice, both interacting self-avoiding walks at the 
collapse temperature and Hamiltonian walks. 

It is possible to perform the loop model to vertex model mapping directly on the 
dilute 0(n)-model to obtain a three-state Temperley-Lieb vertex model (see e.g. Ref. 
|3) with i?-matrix R(u) = PR{u) = 1 + f(u) U where 



U-j — 5i + j : N + i5k+l,N+lQ 2Sl+k 2 , (7.24) 

f(u) = sinh(w)/ sinh(tt + X), N = 3 and q = e x . Integrable boundaries can then be 
derived along the lines of Section 2. While it would be somewhat surprising, it is not 
inconceivable that integrable boundaries other than those in (|7.23 ) could be obtainable. 



8 Discussion 

In this paper we have investigated integrable vertex models and loop models on the 
lattice C. We have applied Sklyanin's theory of K-matrices and transfer matrices 
t(u,uj) to derive integrable weights for various two- and three-state vertex models and 
obtained Bethe ansatz solutions for their eigenvalue spectra. We have also shown how 
some of these integrable vertex models can be interpreted as loop models. It turns out 
that the choice of K- matrices required for this does not usually coincide with the choice 
which leads to quantum group- invariant spin chains. Those 0(n) models amongst the 
integrable loop models have particularly interesting physical applications (in the n —>■ 
limit), being related to various polymer problems in two dimensions. Another set of 
models of interest within the present framework are open boundary versions of the 
restricted solid-on-solid models constructed from loop models [1"3|. These include the 



dilute A-D-E models |43| , |44fl constructed from the 0(n) model of Section 7.1. The 
study of the Bethe ansatz solutions for the loop models found in this paper is currently 
under way and will lead to exact results for various surface critical phenomena. 

The method presented here can, of course, be applied to A^-state models, where 
N > 3, and to related multi-colour loop models ]34]. However, there are still some 



interesting unsolved problems in the case of two-state models. First of all, there is 
the problem of the elliptic eight- vertex model on the lattice C. i^-matrices |J and 
commuting transfer matrices t(u, u)) have been found for this model. Therefore, 
by the method of Section 2, integrable boundary weights can be found. However, 
diagonalization of t(u, u>) is another matter altogether. Because the eight- vertex model 
does not have the property of arrow-conservation across rows, the coordinate Bethe 
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ansatz fails. For the same reason, the algebraic Bethe ansatz presented in Section 3 - 
so successful for the six-vertex model - also fails, because there is no obvious choice 
of pseudovacuum \Q). Baxter's "commuting transfer matrix method" jlj developed 
precisely for the eight-vertex model with periodic boundaries^ might be able to be 
generalized to this open boundary case. 

More or less in the same category is the also interesting problem of diagonalizing the 



transfer matrix t(u } uj) for the six- vertex model with non-diagonal .fT-matrices [31, 36 
The corresponding model on L has the feature of allowing sources or sinks of arrows at 
the boundaries. Presumably some generalization of Baxter's commuting transfer matrix 
method is also required to solve this problem. This is presently under investigation. 

Lastly, but by no means least, integrable vertex models on the lattice C can be 
viewed as convenient regularizations of two dimensional quantum field theories (the so- 
called "light-cone approach" - for a recent review, see Ref. |47| ). In particular, they are 
relevant to the problem of calculating boundary S-matrices for various two dimensional 



field theories with boundaries (see e.g. Refs. |46|, [48], [49], [50|, [51], [53|]). 
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